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Abstract 


The effect of the Gauss-Bonnet term on the existence and dynamical stability 
of thin-shell wormholes as negative tension branes is studied in the arbitrary 
dimensional spherically, planar, and hyperbolically symmetric spacetimes. We 
consider radial perturbations against the shell for the solutions which have the 
Z 2 symmetry and admit the general relativistic limit. It is shown that the 
Gauss-Bonnet term shrinks the parameter region admitting static wormholes. 
The effect of the Gauss-Bonnet term on the stability depends on the space- 
time symmetry. For planar symmetric wormholes, the Gauss-Bonnet term 
does not affect their stability. If the coupling constant is positive but small, 
the Gauss-Bonnet term tends to destabilize spherically symmetric wormholes, 
while it stabilizes hypebolically symmetric wormholes. The Gauss-Bonnet term 
can destabilize hypebolically symmetric wormholes as a non-perturbative effect, 
however, spherically symmetric wormholes cannot be stable. 
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1 Introduction 


Wormhole is a spacetime configuration connecting multiple non-timelike infinities. Worm- 
holes are still hypothetical objects at present but they are certainly realized as solutions 
for the Einstein equations. Theoretical physicists have been attracted to traversable worm- 
holes because they admit an apparently superluminal travel (a spacetime shortcut) as a 
global effect of the spacetime topology PflEIG] and also they are available to make time 
machines HE, !J. (See [I] for a standard text book and [3] for a recent review.) 

As clearly seen in its Penrose diagram, the Schwarzschild vacuum spacetime with positive 
mass certainly possesses the wormhole structure. The static region in this spacetime covered 
by the isotropic coordinates is called the Einstein-Rosen bridge connecting two distinct 
asymptotically flat regions |7J. Of course, the Einstein-Rosen bridge is nothing but a 
static portion of the maximally extended Schwarzschild spacetime, so it represents just an 
instantaneous and non-traversable static wormhole [8j. 

Although the Einstein-Rosen bridge is not satisfactory as a static traversable wormhole, 
such a spacetime can be constructed simply by gluing two Schwarzschild exterior space- 
times pp. There is no singularity in the resulting spacetime and an observer can travel 
from one asymptotically flat region to the other. But unfortunately, this construction re¬ 
quires a massive thin shell at the junction timelike hypersurface which violates the weak 
energy condition. Actually, the requirement of such exotic matters is quite generic for the 
wormhole construction pQ. From a classical viewpoint, exotic matters typically suffer from 
the dynamical instability. However, wormholes might be realized by some quantum effects 
which effectively violate the weak energy condition. For this reason, construction of a dy¬ 
namically stable wormhole in a realistic situation has been a big challenge in gravitation 
physics for a long time. 

In 2011, Kanti, Klcihaus, and Kunz numerically constructed four-dimensional spherically 
symmetric wormhole solutions in Einstein-Gauss-Bonnet-dilaton gravity and showed that 
they are dynamically stable against spherical perturbations |9j. The Gauss-Bonnet term 
non-minimally coupled to a dilaton scalar field appears in the Lagrangian as the ghost-free 
quadratic correction in the low-energy limit of heterotic string theories [lOj. Although this 
Einstein-Gauss-Bonnet-dilaton theory is realized only in ten dimensions, their result gives 
courage and hope toward the construction of wormholes in our universe described by a 
well-motivated effective theory of gravity. Then a natural question arises: Which is the 
main ingredient stabilizing the wormhole, the Gauss-Bonnet term or its dilaton coupling? 

The main purpose of the present paper is to clarify the effect of the Gauss-Bonnet term 
on the dynamical stability. For this purpose, we will study the simplest thin-shell worm- 
hole which is made of its tension |13| . While dynamical stability of thin-shell wormholes 
have been intensively investigated both in general relativity (Einstein gravity) pTJ and in 
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various models of modified gravity [12], this is the best set-up to analyze stability as a pure 
gravitational effect because such a thin shell does not suffer from the matter instability. 


In Einstein gravity, such thin-shell wormholes have been fully investigated by two of the 
present authors in the spacetimes with a maximally symmetric base manifold in arbitrary 
dimensions, namely arbitrary-dimensional spacetimes with spherical, planar, or hyperbolic 
symmetry [EE]. In the vacuum case, such thin-shell wormholes are stable against radial 
perturbations only in the hyperbolically symmetric case with negative mass in the bulk 
spacetime 0- 

In the present paper, we will study the same system with the Gauss-Bonnet term but 
without a dilaton in the Lagrangian. Since the Gauss-Bonnet term becomes total derivative 
and does not affect the field equations in four or less dimensions in the absence of the non- 
minimal coupling to a dilaton, we will consider five or higher-dimensional spacetimes. In 
comparison with the general relativistic case, the equation of motion for the shell is much 
more complicated. For this reason, although thin-shell wormholes have been investigated 
in Einstein-Gauss-Bonnet gravity by many authors [15], the stability analysis has not been 
completed yet even against radial perturbations. 

The outline of the present paper is as follows. In section II, we derive the equation of 
motion for a shell and the basic properties of the static shell are reviewed. In section III, 
we study the existence and stability of static thin-shell wormholes in the cases except for 
k — — 1 with m < 0. Section IV is devoted to performing the pictorial analysis to study the 
same problem in the case of k — — 1 with m < 0. Our results are summarized in section V. 
A detailed derivation of the equation of motion for a thin shell is summarized in appendix 
A. Our basic notation follows 0- The convention for the Riemann curvature tensor is 
[V p , Vo-] V p = R^vpaV 11 and = R p p, pu . The signature of the Minkowski metric is taken 
as diag(—, +,+,••• ,+,+), and Greek indices run over all spacetime indices. We adopt the 
units in which only the d-dimensional gravitational constant Gd is retained. 


2 Preliminaries 

2.1 Einstein-Gauss-Bonnet gravity 

In the present paper, we consider d(> 5)-dimensional Einstein-Gauss-Bonnet gravity in 
vacuum, of which action is given by 



( 2 . 1 ) 
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where Hd := \Z&nGd and A is the cosmological constant. The Gauss-Bonnet term L GB is 
defined by the following special combination of the Ricci scalar R. the Ricci tensor R IW and 
the Riemann tensor R^upa- 

L gb := R 2 - 4 R^RT + R^RT**. (2.2) 

The Gauss-Bonnet term appears in the action as the ghost-free quadratic curvature correc¬ 
tion term in the low-energy limit of heterotic superstring theory in ten dimensions (together 
with a dilaton) p[UJ. In this context, the coupling constant a is regarded as the inverse 
string tension and positive definite. For this reason, we assume a > 0 throughout this 
paper. In addition, we put another conservative assumption 1 + 4bA > 0, where 

2A 

A:= (d-l) ( rf-2) ’ (<i-3)(d-4)a, (2.3) 

so that the theory admits maximally symmetric vacua, namely Minkowski, de Sitter (dS), 
or anti-de Sitter (AdS) vacuum solutions. Although there exists a maximally symmetric 
vacuum for 1 + 4bA = 0, we don’t consider this case for simplicity. 


Variation of the action (2.1) with respect to the metric g A 
Einstein-Gauss-Bonnet equations: 


gives the following vacuum 


G\ + oiH\ + A W = 0, 


(2.4) 


where 

Gfiv • (2.5) 

Hpu :=2 (rR^v - 2 Rp Q R a u - 2 R ad R^p + R,, nB l Ruap-^j - ^g^L GB . (2.6) 

The tensor H /w obtained from the Gauss-Bonnet term does not give any higher-derivative 
term and Hp U = 0 holds for d < 4. As a result, Einstein-Gauss-Bonnet gravity is a second- 
order quasi-linear theory as Einstein gravity is. 


2.2 Bulk solution 

We will study the properties of thin-shcll wormholes in Einstein-Gauss-Bonnet gravity. 
Such wormhole solutions are constructed by gluing two bulk solutions at a timelike hyper¬ 
surface. 

In the present paper, we consider the d-dimensional vacuum solution with a maximally 
symmetric base manifold na as the bulk solution, of which metric is given by 

= g^dx^dx" = - f(r)dt 2 + /(r) _1 dr 2 + r 2 'y AB dz A dz B , (2.7) 
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where z " 4 and 7 ab (A, B — 2, 3, • • • , d — 1) are the coordinates and the unit metric on the 
base manifold and 

(W i+ ^ + *4 <2 - 8) 

Here k = 1,0, — 1 is the curvature of the maximally symmetric base manifold corresponding 
to the spherical, planar, and hyperbolically symmetric spacetime, respectively, m is called 
the mass parameter. 


The expression of the metric function (2.8) shows that there are two branches of solutions 
corresponding to the different signs in front of the square root. The branch with the minus 
sign, called the GR branch, allows the general relativistic limit a —> 0 as 


f(r) = k 


m 

r d— 3 


A r 2 


(2.9) 


On the other hand, the metric in the branch with the plus sign, called the non-GR branch, 
diverges in this limit. In the following section, we will consider the bulk solution only in 
the GR branch as a conservative choice. 


The global structure of the spacetime (2.7) depending on the parameters has been com¬ 
pletely classified 


pl]. There are two classes of curvature singularities in the spacetime. 
One is the central curvature singularity at r = 0. Since we assume a > 0 and 1 + 4<5A > 0, 
the interior of the square root becomes zero at some r = r^(> 0) for negative m. This 
corresponds to another curvature singularity called the branch singularity and the metric 
becomes complex at r < 77 . In this case, the domain of the coordinate r is r G ( 77 , 00 ). 


The spacetime has a Killing horizon at r = 77 satisfying /(rh) = 0 depending on the 
parameters. I 11 order to construct static thin-shell wormholes, the bulk spacetime needs to 
be static. For this reason, we consider the bulk solution ( |2.7 ) in the domain r G ( 77 , 00 ) 
if there is no branch singularity and r G (max{7y,, rh}, 00 ) if there is a branch singularity. 
We define the future (past) direction by increasing (decreasing) direction of t. 


2.3 Equation of motion for a thin shell 


A thin-shell wormhole spacetime is constructed by gluing two bulk spacetimes (2.7) at a 


timelike hypersurface r = a. Here the bulk spacetimes are defined in the domain r > a{> 
77 ) and may have different parameters. Then the junction conditions, which are the field 


equations (2.4) in the distributional sense, tell us the matter content on the thin shell at 
r = a. Finally, the equation of motion for the shell is obtained as a closed system when an 
equation of state for the matter is assumed. 
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The junction condition in Einstein-Gauss-Bonnet gravity is given by 


[tfj]± - 6‘jlKU + 2a(3[jy ± - <q[J] ± - 2 P‘ kjl {K M } ± ) = -K 2 d S‘ p (2.10) 


where i,j — 0,1, - - - , (d — 1) are indices for the coordinates on the timelike shell |[[9| 20] , 
Here we have introduced 


[A] ± := A+ - A'" , 


( 2 . 11 ) 


where X ± is the quantity A" evaluated either on the + or — side of the shell. In the junction 
conditions (2.10), K‘j is the extrinsic curvature of the shell and K := h^Kij, where h x j is 
the induced metric on the shell. Other geometrical quantities are defined by 


Jij ■= ^{2KK ik K k j + K kl K kl K ij -2K ik K kl K lj -K 2 K ij ) , (2.12) 

Pikji • T^ikji + 2h i{l K j]k T 2h k \ J 7Zi} l T TZ1 i t /p j/,. , (2.13) 


where lZ X jki , 7and 7Z are the Riemann tensor, Ricci tensor, and Ricci scalar on the shell. 
Pijki is the divergence-free part of the Riemann tensor 7 Zijki satisfying DiP l jki = 0, where 
D{ is the covariant derivative on the shell. Lastly, S l j is the energy-momentum tensor on 
the shell, which satisfies the conservation equations D l S l J = 0. 


A static thin-shell wormhole is realized as a static solution for the equation of motion. 
However in general, a is not constant but changes in time, representing a moving shell. In 
such a case, a may be written as a function of the proper time r on the shell as a = a(r). 


Now let us derive the equation of motion for the shell. We describe the position of the 
thin shell as r = a(r) and t — T(t) in the spacetime (2.7) and assume the form of S l j as 

S l j = diag (~p,p,p, ■ ■ ■ ,p) + diag(—cr, -a, -a, • • • , -a,). (2.14) 

This assumption means that the matter on the shell consists of a perfect fluid and the 
constant tension cr of the shell, where p and p are the energy density and pressure of 
the perfect fluid. Assuming the Z 2 symmetry for the bulk spacetime, we write down the 
junction conditions (2.10) as 

1 


7^d(P + T> = 


(d — 2)fT f , 2d 


V 


.a 2 


3 k f 


1 I 2^7 H , 77 


(2.15) 


-^ 2 d (jp-cj) = -^p \7 + T7 + ( d-3 )f^ + i 


fT {a 2a 
2aa f d — 5 / a 2 


fT V 


/ 


3 \a 2 


a- 


d 2 


3 k f 




a z 


a 2 k f \ a. f (k f 

2— + — -I -^) _ + X _ (~-2 

a z cr a z j a 2a V a z a z 


(2.16) 
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where / = /(a). A dot and a prime denote the derivative with respect to r and a, 
respectively. (See Appendix |A] for the details of derivation.) The above equations give the 
equation of motion for a thin shell in Einstein-Gauss-Bonnet gravity. 


In order to obtain the equation of motion in a closed system, an equation of state for 
the perfect fluid is required. One possibility is the following linear equation of state with a 
constant 7 : 

p=( 7 -l)p. (2.17) 


With this equation state, the energy-conservation equation on the shell 
as 

P — ~(d— 1 )(p + p)~, 
a 

is integrated to give 

_ Po 
P a 0 _ 1 ) 7 ’ 


where po is a constant. 


0 , written 
(2.18) 

(2.19) 


2.4 Effective potential for the shell 


The dynamics of the shell governed by Eqs. (2.15) and (2.16) with an equation of state (2.17) 


can be discussed as a one-dimensional potential problem. Then the shape of the effective 
potential V(a) for the shell determines the stability of static configurations, namely the 
static wormholes. 


Let us derive the effective potential V(a). Squaring Eq. (2.15) and using Eq. (A.24), we 
obtain 


/ a 2 \ f 2 / d 2 3k f 


^ L 2 + a 2 Ji 1 + 3°Va 2 + a 2 a 2 


( 2 . 20 ) 


where 


G(a ) 2 : 


2 «d(p( a ) + °) 


4(d — 2 ) 2 ' 


( 2 . 21 ) 


This is a cubic equation for a 2 . The position of the throat a = a 0 for a static wormhole is 
obtained by solving the following algebraic equation for ao: 


„(3k _fo 

0 a o l 3 a U 4 


( 2 . 22 ) 


where /o := /(ao) and Gq := ^( a o) 2 - 









For convenience, we define 


,, , 4am 

A{r) :=1 H- 3 —r + 4aA 


~d— 1 


(2.23) 


with which the metric function (2.8) in the GR branch is written as 

„2 


f(r) =k+ ^{ 1 ~ V^r)) • 


(2.24) 


A > 0 is required for the real metric and the absence of branch singularity. In the GR 
branch, because of the exitence of the squre root in Eq. (2.24), the following inequality 
holds: 


r 2 + 2 ak — 2 af(r) > 0 , 


(2.25) 


which will be used later. 


Actually, Eq. (2.20) admits only a single real solution for a 2 : 

a 2 = —V(a), 


(2.26) 


which has the form of the one-dimensional potential problem. The effective potential V (a) 
is defined by 


V(a) := f(a ) — J(a)a 2 , 


where J(a) is defined by 

(B(a) - A(ay/°)'- 


J(a) 


4 aB(a) 


One can see B > A 1//2 . Q 2 can be expressed in terms of A and B 


as 


0 2 = 


(B 3 - A 3 / 2 ) 2 
36 aB 3 


(2.27) 


(2.28) 


B(a) :=|l 8 dO(a ) 2 + A(a) 3/2 + 6 y/afl(a) 2 (9af2(a ) 2 + A(a) 3 / 2 )| . (2.29) 


(2.30) 


2.5 Existence conditions for static shell 


Here we summarize the existence conditions for a static shell located at a = oq. Equa¬ 


tion (2.26) is interpreted as the conservation law of mechanical energy for the shell. By 
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differentiating Eq. (2.26) with respect to r, we obtain the equation of motion for the shell 


as 


a = --V'(a) 


(2.31) 


From Eqs. (2.26) and (2.31), a 0 is determined algebraically by U(a 0 ) = 0 and V'(a 0 ) = 0. 


In addition, a 0 must satisfy A(a 0 ) > 0 and f(a 0 ) > 0. The latter condition f(a 0 ) > 0 is 
called the horizon-avoidance condition in Ref. m, which simply means that the position 
of the throat is located in the static region of the spacetime. Actually, this condition is 


always satisfied because we have U(a 0 ) = 0 and Eq. (2.27) implies /(a) > V(a). 


2.6 Negative energy density of the shell 


In closing this section, we show that the energy density on the shell p + cr must be negative 


for static wormholes. The condition p + a > 0 and Eq. (2.15) with a = a 0 (> 0) yields 

4 ak 


«o< 


(2.32) 


2 + a/A) 

where A 0 := A(a 0 ). Clearly, this is not satisfied for k = 1,0 under the assumption a > 0. 


For k = —1. Eq. (2.32) gives 


«o< 


4 a 


2 + a/A) 


< 2 a 


(2.33) 


and this is not satisfied because there is a constraint a q > 2 a for the throat radius in 


the GR branch, which can be shown from the combination of Eq. (2.25) and /(ao) > 0. 
Now we have shown that the energy density on the shell is negative in the physical set up 
considered in the present paper. 


3 Static thin-shell wormholes made of pure tension 


In the present paper, we analyze stability of the static shell located at a = against radial 
linear perturbations. The shape of the effective potential V ( a ) determines the dynamical 
stability of the static shell, as explained below. By Eqs. (2.26) and (2.31), a 0 satisfies 


V(ao) = 0 and V'(ao) = 0. Using U(ao) = V\aa) = 0, we obtain the Taylor expansion of 
the potential V (a) around a = ao as 


V (a) — -V"(cio)(a — a o ) 2 + 0((a — ao) 3 ). 


(3.1) 
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The stability condition against radial perturbations for the shell is then given by 

V"{a Q ) > 0 . (3.2) 

Hereafter we will consider the case without a perfect fluid on the shell (p = p = 0) and 
assume a < 0. The resulting static thin-shcll wormholes are made of the pure (negative) 
tension a and satisfy the null energy condition. This simplest set up is preferred by the 
minimal violation of the energy conditions. A technical advantage in this set up is the 
constancy of fl 2 . 


3.1 Existence conditions 


The location of the static wormhole throat a = ao is determined by the following algebraic 


equation obtained by eliminating a from Eqs. (2.15) and (2.16): 


h-hV 1 + 2a k ~ fa 

a 0 2 


Aakf 0 


= 0, 


(3.3) 


where f' Q := /'(a 0 ). In the limit of a —>• 0, Eq. (3.3) reduces to the corresponding equation 
in Einstein gravity [3]. The explicit form of Eq. (3.3) is 


kal\^A 0 = 2kai — --b 4 ak 2 . 

0 V 0 2 at 5 


(3.4) 


This equation shows that m = 0 is required for k = 0 and then ao is totally undetermined in 


the domain where both A 0 > 0 and / 0 > 0 hold. The metric function (2.8) with m — k — 0 


shows that the horizon avoidance condition f 0 > 0 is satisfied only for A < 0. Stability of 


this wormhole will be investigated in Section 3.5.1 


On the other hand, for k = 1 (—1), the left-hand side of Eq. (3.4) is positive (negative) 
and hence the throat radius must satisfy 


kl 2 ka. 


2 (d — 1 )m 


2a! 


d— 5 


+ 4a > > 0 




For k 7 ^ 0, squaring Eq. (3.4) gives the following algebraic equation for a 0 : 


(3 — 4dA)ao^ 3 ' ) + 16aka^ d + 16a 2 a^ d 5 ' ) 

— 2 {d — l)kmaQ~ 3 


Adama q 5 + -{d — l) 2 m 2 = 0. (3.5) 


Static wormhole solutions with the throat radius ao must satisfy Eq. (3.5) and also several 
constraints. 
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The first constraint is the inequality (A). Another constraint comes from Eqs. (3.34) 
and (3.35). Eliminating /q, we obtain 


A =i{^- 2 « +2a h- 


(3.6) 


Since the right-hand side of Eq. (3.6) must be positive, we have a necessary condition for 
physical solutions: 


k{ ^ J 1 } m - 2 (kal + 2d))> 0. 


^ d —5 




(B) 


We do not have to impose the condition A 0 > 0 to avoid the complex metric in the 


bulk spa cetime because any real solution of Eq. (3.5) satisfies it. This is shown as follows. 
Eq. (3.5) is solved for m/a g~ 5 as 


2 \ 2k(d — l)ag + 4 da ± \j(d — 1) 2 (1 + 4aA)ag + 16 k(d — l)aag + 16(2d — l)a 2 


^d —5 


(d-iy 


Substituting this into A 0 = 1 + 4dA + 4<3m/ao 1 , we obtain 

U ± 8 a\JU — 16b 2 


An 


(d — l) 2 ao 


(3.7) 


(3.8) 


where 


U :=(d — 1) 2 (1 + 4«A)aQ + 16a { k(d — l)ag + 2 da} . (3.9) 

For any real solution of Eq. (3.5), the interior of the square root in Eq. ( |3.7 ) is non-negative, 
which gives the following lower bound of A: 

(d — l) 2 ag + 16 (d — l)aka% + 16(2 d — l)b 2 


A > - 


4 (d, — l) 2 aag 


(3.10) 


This inequality implies U > 16d 2 and hence U is positive. Therefore Aq with the plus sign 


in Eq. (3.8) is positive. Positivity of Aq with the minus sign is shown by direct computations 


without using the inequality (3.10). 


3.2 Non-existence for k = 1 with m < 0 and k = —1 with m > 0 


It is shown that there is no static thin-shcll wormhole for k — 1 with m < 0 and k = — 1 


with m> 0. For the proof, we use Eq. (2.26) in the following form: 


d In a 
dr 


+ V(a) — 0, 


(3,11) 
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where 


V(a) := ^ - J (a), 
cr 

There is no static solution if V(a) is monotonic. V' is computed to give 


V' 


1 [ 8k ®B 2 . (B 2 

4 aB 2 \ a 3 1 


+ BA' 


The following expressions; 


.. . 4 (d — 1 )am 

(«) =-^, 

. ^(a) 1/2 A'(a) / _ 3afi 2 _ 

2 B 2 f \/“^ 2 (9a!2 2 + /1(a) 3 / 2 ) 


(3.12) 


(3.13) 


(3.14) 

(3.15) 


imply B' < 0 (> 0) and A' < 0 (> 0) for m > 0 (m < 0) with equality holding for m = 0. 
Together with the facts A, B > 0 and B > A 1 / 2 , it is shown that V' is negative definite for 
k — 1 with m < 0 and positive definite for k = — 1 with m > 0. 


3.3 Stability criterion 

As explained at the beginning of this section, stability of a static thin-shell wormhole is 
determined by the sign of V"(ao). In this subsection, we will derive V"(ao) in a convenient 
form to prove the (in)stability. 


3.3.1 General relativity 


First let us consider Einstein gravity as a simple lesson. In the general relativistic limit 


a —> 0, Eq. (2.20) reduces to 


IT 


f( a 


a 2 + a 2 ' 


(3.16) 


By solving Eq. (3.16) for a 2 , we define a potential V(a) of the conservation law of the 
one-dimensional potential problem. Then we directly calculate the second derivative of 
V(a). However, without such direct calculations, in principle we can derive the form of 
V"(ao) by operating a systematic method below, which can be applied also in more general 
theories of gravity. 
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Suppose we get a master equation as a 2 + V(a) = 0. By this master equation, a 2 in 
Eq. (3.16) is replaced by —V(a) to give 


^2 = f{ a ) -V(a) 

a 2 

Differentiating this equation twice, we obtain 

0 =«(/' - V') - 2(/ - V), 
0 =a (/" - V") - (}’ - V'). 


(3.17) 


(3.18) 

(3.19) 


In Einstein gravity, the metric function is 

/(«) 

which satisfies 

(d - 3)(/c - /) - A (d - 1 )a 2 

J (a) = -, 

a 

A (d - 1)0 - 4)a 2 - (fc - f){d- 2)(d - 3) 


(3.20) 

(3.21) 

(3.22) 


Substituting Eq. (3.21) into Eq. (3.18) and evaluating it at a 
E'Oo) = 0, we obtain 


ao satisfying V(ao) = 


/o(:= /(ao)) = ~ 


(3.23) 


Combining this with Eq. (3.20), we obtain the algebraic equation to determine ao'- 

2k 


m 


d- 1 


,< 2-3 ‘ 


(3.24) 


For k — 0, Eq. (3.24) requires m = 0 and ao is totally undetermined. For k = 1(—1), 
Eq. (3.24) requires m > (<)0 and the throat radius ao is given by 

(d — l)m\ 1 ^ d ~ 3 ' > 


ao — 


2k 


(3.25) 


As seen in Eq. (3.25), A does not contribute to the size of the wormhole throat. However, 
it appears in the horizon-avoidance condition f 0 > 0. Equation (3.23) shows that f 0 > 0 
requires A < 0 in the case of k = 0, —1. In the case of A = 0, f 0 > 0 is satisfied only for 
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k — 1. In the case of A > 0 and k — 1, /o > 0 gives a constraint ao < cic ; on the size of 
the wormhole throat, where 


(GR) ( (d 3)fc \ / 

c ' V(d-i)A; 

On the other hand, in the case of A < 0 and k — — 1, /o > 0 gives ao > a. 


(GR) 


(3.26) 


Combining 


this inequality with Eq. (3.25), we obtain the range of the mass parameter admitting static 


/ /~i id \ f’C'R'l 

wormhole solutions; 0 < m < ; for k = 1 with A > 0 and m < ml {< 0) for k = —1 

with A < 0, where 


m, 


(GR) 


2k ((d - 3)k\ {d ~ 3)/2 


d-l\(d-l)A 


) 


(3.27) 


In Einstein gravity, a simple criterion for the stability of static solutions is available. 
Substituting Eqs. (3.21) and (3.22) into Eq. (|3.19[), evaluating them at a = ao, we obtain 


(d-l)(d-3)m 

^ («o) =- -jzi - 

a 0 

2(d — 3)k 

2 ) 
a 5 


(3.28) 


where we used Eqs. (3.20) and (3.24). This simple expression clearly shows that the worm- 

Existence and stability of static thin-shell 


hole is stable only for k = — 1 with m < 0 
wormholes in Einstein gravity are summarized in Table [lj 


Table 1: The existence and stability of Z 2 symmetric static thin-shell wormholes m ade o f 
pure negative tension in Einstein gravity, where ac GR ' ) and mf R ^ are defined by Eqs. (3.26) 
and (3.27), respectively. 




Existence 

Possible range of a 0 

Stability 

k = 1 

A > 0 

n . 7 (GR) 

0 < m < ml 

rTZ 7 (GR) 

0 < a 0 < a c 

Unstable 

A < 0 

m > 0 

ao > 0 

Unstable 

k = 0 

A > 0 

None 

- 

- 

A < 0 

rri — 0 

ao > 0 

Marginally stable 

k — —1 

A > 0 

None 

- 

- 

A < 0 

m < mf K, (< 0) 

a 0 > ac 

Stable 
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3.3.2 Einstein-Gauss-Bonnet gravity 


Although it is more complicated, we can play this game in Einstein-Gauss-Bonnet gravity 
in a similar manner. Replacing a 2 by —V(a) in the master equation (2.20), we obtain 


ft 2 


/(a) - V(a) f 2a(—2V (a) + 3k - /(a)) 


1 + 


3a 2 


(3.29) 


In Einstein-Gauss-Bonnet gravity, the metric function is 

-.2 


which satisfies 


/'(«) 

/"(«) 




(<d — 5 )a(k — f) 2 + (d — 3)a 2 {k — /) — A (d — l)a 4 


a{a 2 + 2 a(k — /)} 


L{a) 


a 2 {a 2 + 2a (k — /)} 3 ’ 


where 


(3.30) 


(3.31) 

(3.32) 


L(a) :—2(d — l) 2 aA 2 a 8 — Aa 4 (d — l)|l2a 2 (fc — /) 2 + 12 aa 2 {k — /) — (d — 4)a 4 
-(/; — /){ 2(d -3)(d- 5)a 3 (k - f ) 3 + 4(d 2 - 8d + 13)dV(fc - /) 2 


+ 3(d — 2){d — 5)aa 4 (k — f) + {d — 2){d — 3 )a e 


(3.33) 


Equation (3.30) gives 


m = a d 3< ( —A a 2 + (k — f(a )) + aa 2 {k — f(a )) 


(3.34) 


Differentiating (3.29) and evaluating at a = a 0 , we obtain 

2 {(d — l)a 2 + 2fc(d + l)d}/o + (d — l)Aag — k{(d — 3)oq + (d — 5)dfc} 


/« 


o 


(d — l)a 


(3.35) 


where we used Eq. (3.31). This equation reduces to Eq. (3.23) for a —» 0. Equation (3.35) 


will be used to replace /q (p = 2,3,4, • ■ ■) by f 0 . 
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Differentiating Eq. (3.29) twice and using Eqs. (3.31) and (3.32), we finally obtain E"(a 0 ) 
in a rather compact form: 


V"(a 0 ) = 


2kP(ao ) 


a o( a o + 2fca + 2a/o)(ao + 2 ka — 2d/ 0 ) ’ 


P(a 0 ) := 4d 2 / 0 < 6k — (d — 3)/ 0 f+(aQ + 2fcb)< (d — 3)ap + 2(d — 5)/cb >, 


(3.36) 

(3.37) 


where we have eliminated A by using Eq. (3.35). This expression reduces to Eq. (3.28) for 


a —y 0. Because of Eq. (2.25), the denominator in the expression of V"(ao) is positive and 


therefore the sign of the function P(ao) determines the stability of the shell. 


3.4 Effect of the Gauss-Bonnet term on the stability for a/a 2 E <C 1 


Before moving onto the full-order analysis, let us clarify how the Gauss-Bonnet term affects 
the stability of thin-shell wormholes in the situation where a is small. 


In the general relativistic limit a —>• 0 with k = ±1, Eq. (3.4) gives 


d -3 _ (d — l)mk __ d _ 


a o ~ 


-• OjT 


(3.38) 


This is the static solution in Einstein gravity which requires mk > 0 [Ti] . Now we obtain 
the static solution for e := a/a? E 1 in a perturbative method. We expand ao in a power 
series of e : 


&0 — OE + ®(l) e T" ®(2) e + • • 


(3.39) 


Substituting this expression into Eq. (|3.4|) and expanding it in a series of e, we obtain 

(3.40) 


2A(d — l)a| — 4(d — 2)k 

a (i) —-- a E- 


(d — l)(d — 3) 

This allows us to derive the expansion of Eq. ( |3.36[ ): 

T//// , ^ 2(d — 3)k t 4(d — 3)fca(i)e : 8 k 2 e 

* (^0) - O T o o 

2(d — 3)k 8ke (d — 3 , ~ 2 

' --I j—^k ~ Aag 


(3.41) 


The first term of Eq. (3.41) coincides with Eq. (3.28) and inside the bracket of the 


second term is positive because of Eq. (3.23). Hence we arrive a simple conclusion about 


the effect of the Gauss-Bonnet term for small b; it destabilizes wormholes in the spherically 
symmetric case (k = 1), while it stabilizes in the hyperbolically symmetric case (k = —1). 
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3.5 Stability for k = 0,1 


In this subsection, we will prove (in)stability of the static thin-shell wormhole in the frame¬ 
work of Einstein-Gauss-Bonnet gravity. We are going to study the sign of V"(ao) given by 
Eq. (3.36) for k = 0 and k — 1. Because the analysis is much more complicated in the case 
of k — —1, it will be treated in the next section. 


3.5.1 k — 0 with m = 0: Marginally stable 


The analysis for k = 0 is very simple. For k = 0, Eq. (3.4) 
undetermined. Therefore, any size of the static wormhole 
This is consistent with the fact that Eq. (3.36) gives V"(ao ) 
marginally stable. 


gives m = 0 and ao is totally 
throat is allowed for k = 0. 
= 0, namely the wormhole is 


3.5.2 k = 1 with m > 0: Unstable 


In Section 3.2, we have shown that there is no static wormhole for k — 1 with m < 0. In 


the present paper, we don’t clarify the parameter region with positive m admitting static 
wormhole solutions because they are all dynamically unstable in any case. In Fig. [lj we 
plot the profile of V(a ) with k = 1 and m > 0, in which there is a local maximum. This 
implies that the corresponding static solution is unstable. We will prove this analytically. 


For k = 1, positivity of P(ao) in Eq. (3.36) means instability of the static wormhole. 


Using the inequality (2.25), we evaluate the lower bound of P(a 0 ) as 


P(a 0 ) =4a 2 / 0 ^ 6 — (cl — 3)/ 0 J> + (ao + 2a)|(d — 3 )a^ + 2 (d — 5 )a 


>4cr/ 0 < 6 — (d — 3)/ 0 >+2a/ 0 < (d — 3)aQ + 2 (d — 5)b 


—2 a fo{ 2 (d — 3 )a ( ^ — ) +2(d + l)a 


>2 afo{ -2(d - 3)a + 2 (d + l)d \= 16 a 2 f 0 > 0. 


(3.42) 


Therefore, the wormhole is dynamically unstable. 
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V(a) 



Figure 1: The potential V(a) for d = 5,6,7 in Einstein and Einstein-Gauss-Bonnet (EGB) 
gravity with k = 1 , a — 0 . 02 , m — 1, A = 1 and cr = — 0 . 1 . 


4 Stability for k = — 1 with m < 0 


In the present section, we will provide the stability analysis for k = —1. 
shown that there is no static wormhole for k = — 1 with m > 0 in Section 322 
discuss the case with m < 0. 


Since we have 

we will 


4.1 Preliminaries for pictorial analysis 


For our purpose, we introduce x := cIq and y := m/a q 5 , with which Eq. 
as h(x, y ) = 0, where 

h(x, y) :=(3 — AaA)x 2 — 2 (d — 1 )kxy + — l) 2 y 2 + 16 akx — Aday + 16d 2 . (4.1) 

We adopt a pictorial analysis in the (x,y) plane in the domain of x > 0 and y < 0. 

For 3 — 4bA ^ 0, h(x, y) — 0 is solved to give x = x±(y), where 


(3.5) is written 


x±(y) := 


2 k{(d - 1 )y - 85} ± y/Z(y) 
2(3 - 4bA) 


(4.2) 
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The function Z(y) in the above expression is defined by 


Z(y) :=4|(d — 1 )y — 85 > —(3 — 4bA)<j(d — 1 ) 2 y 2 — 16 day + 64b 2 

—(d — 1) 2 (1 + 4bA)y 2 — 16b(4dbA + d — 4)y + 64b 2 (l + 4bA). (4.3) 

In the limit a —> 0, we obtain 


lim x+(y) 

a —>0 


k(d-l)y 


lim x _ (y) = 

ct —>-0 


k(d~ l)y 
6 


(4.4) 


Among these two, only the former satisfies Eq. (3.4) with a = 0. For this reason, we will 
focus only on x + (y) hereafter because X-(y) does not admit the general relativistic limit. 


4.1.1 Physical domain of solutions 


Solutions of Eq. (3.5) are realized as intersections of h(x,y) = 0 with y = m/x^ d ~ 5 ^ 2 in 
the (x, y) plane. In addition, they must be located in the physical domain where all the 
constraints on the solutions are satisfied. 


The inequality (^4) gives a constraint between x and y for physical solutions: 

(d-l)k 


x > 


-y - 2ak =: x min (y). 


Also, the inequality (|Sj) gives another constraint: 

(d-l)k 


—2 ka < x < 


y - 2 ka =: x max (y), 


(4.5) 


(4.6) 


where the left inequality comes from a q + 2 ka > 0. In summary, physical solutions must 
be located in the domain of x m ; n (y) < x < x mgbX (y) and y < 0. Since a static solution is 
realized as an intersection of the hyperbola with y = m/x( d ~ 5 ^ 2 in this domain, the number 
of static solutions depend on the value of m. (See Fig. [2] as an example.) 


4.1.2 Stable domain of solutions 


Substituting Eq. (3.6) into Eq. (3.37), we write P(a 0 ) as a function of x and y as 


P(x, y) =k< (d — 1 )(d — 3 )y — 16b >x 


— (d — 3 ){d — 1 ) 2 y~ — 8 d(d — 1 )ay + 128a^ 


(4.7) 
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X 



Figure 2: The (x,y) plane for —1 < 4bA < —( 2d — 5)/(2d — 1) with k = —1, a > 0, 
and d — 6. The thick hyperbola consists of x — x + (y) and x = X-(y ), while dashed lines 
consist of x = x max (y) and x = x m i n (y). Thin curves correspond to y = m/x with 
three different values of negative m. It moves to the left as m(< 0) decreases. Since a 
static solution is realized as an intersection of the hyperbola with y = m/x( d ~ 5 '>/ 2 in the 
shadowed region, the number of static solutions depend on the value of m. 
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The stable domain of solutions are given by P(x,y ) > 0 in the (x,y) plane. The curve 
P(x, y) = 0 representing marginal stability is given by 


(d — 3 )(d — 1 ) 2 y 2 — 8 d(d — 1 )ay + 128a 2 
Ak{(d — 1 )(d — 3 )y — 16d} 


= : x P (y). 


We have xp(0) = —2 ka and 


lim xp (y) 

y —>—oo 


d — 1 



2 da 

k(d — 3) 


(4.8) 


(4.9) 


Since P(x, y) is an increasing function of x in the negative domain 
P(x min (y),y) =2(d- 1 )ay < 0, 

P{x^ 3 »,y) = ^ (<*-!)»{»«-(<<-!)(«<-%} > 

xp( y) satisfies x min (y) < x P (y) < x mSiX (y). (See Fig. 0 as an example.) 

Now we are ready to perform the stability analysis. We will treat 
and — 1 < 4bA < 0, separately. 


of y and we have 

(4.10) 

0, (4.11) 

the cases of A > 0 


4.2 Non-existence for A > 0 

In this subsection, we treat the case of A > 0. Similar to the general relativistic case, 
Einstein-Gauss-Bonnet gravity also does not admit static thin-shcll wormholes in this pa¬ 
rameter region. 


For 4aA > 3, Eq. (4.2) shows that x+ < 0 < x_ holds in the domain of y < 0 with 


k = — 1. Since x = x+(y) does not satisfy the necessary condition x > 0, there is no static 
solution in this case. 


In the special case of 4aA = 3, h(x, y) = 0 is solved to give 

. (d — 1 ) 2 y 2 — 16 day + 64b 2 , . 

x (y) = ——(> °)» 


which satisfies 


8 k{(d — 1 )y — 8a} 


*-.(») - <v) = — > o, 


8 k{(d — 1 )y — 8a} 


(4.12) 


(4.13) 


where x mi n (y) is defined in Eq. (4.5). Because the necessary condition x > .T min is not 
satisfied, there is no static solution for 4bA = 3. 
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Lastly, in order to show the non-existence for 0 < 4bA < 3, we will use the following 
fact: The sign of x max (y) — x + (y) is definite in some domain of y if the curves x + (y) and 
i ma x (y) do not intersect and x + (y) is continuous there. Actually, x +(y) is continuous in 
the negative domain of y for 4bA > — ( 2d — 5)/(2d — 1) because Eq. (4.3) shows that Z(y ) 
is non-negative then. 

It is also possible to show that x + (y ) and x m & x (y) do not intersect in the negative domain 
of y. For A = 0, x + (y ) = x max (y) is solved to give 

»=^s- (414) 

which is positive. For 0 < 4bA < 3, the solution is 


4 (d — l)aA + d — 3 ± 


V 


l)aA + d — 3 > — A(d — l) 2 aA.(l + 4bA) 


(d- 1) 2 A 


A(d — l)aA + d — 3 ± \ A(d — 1 )(d — 5)aA + (d — 3 ) 2 


• 2/c(±)j 


(4.15) 


(d - 1) 2 A 

where inside the square-root is positive for 4bA > —1. By direct calculations, both y c {+) 
and y c (~) are shown to be positive for A > 0. (We note that y c (+) < 0 and y c {~) > 0 are 
satisfied for A < 0.) 


We have shown that the sign of x + (y ) — x max (y) is definite in the domain of negative y 
for A > 0. This sign is actually negative, as shown below. From the following expression; 

k(d — 1)(1 — AaA)y + Aka(l + 45A.) — 4 J Z(y) 


-^max (: v ) -x+(y) 


2(3 - 4bA) 

-Aka{(d — 1 )y — 4<3}A + k(d — 1 )y + Aka — 4 / Z(y) 


we obtain 


*^max (0) -1+(0) 


2(3 - 45A) 

2ka(l + 4dA) — Aoi\J 1 + 4dA 


(4.16) 


(4.17) 


3 -4dA 

which is negative for —1 < 4bA < 3. Because x + (y) — x max (y) is continuous, it is concluded 
that x + (y) > x mSLX (y ) is satisfied in the negative domain of y for 0 < 4bA < 3. 


4.3 Pictorial analysis for —1 < 4<SA < 0 

Now we focus on the case of —1 < 4bA < 0. In this case, there exit static solutions but 
their existence and stability depend on the parameters in a complicated manner. We will 
clarify them by a pictorial analysis. 
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4.3.1 Geometric shape of h(x,y) = 0 


In the case of —1 < 4bA < 0, h(x,y) = 0, where h(x,y) is defined by Eq. (4.1), is a 
hyperbola in the (x, y) plane in general. In order to understand the shape of the hyperbola, 
we present 


h( 0, y) =-(d — 1 ) 2 y 2 — Aday + 16b 2 , 
h(x, 0) =(3 — 4bA)a; 2 + 16 akx + 16b 2 . 
h(0,y) = 0 has two positive solutions y = y\(> 0) and y = t/ 2 (> Vi), where 


V i := 


8 a(d — y/2 d — 1) 


V2 ■ = 


8 a(d + \j2d — 1) 


(4.18) 

(4.19) 

(4.20) 

(4.21) 

(4.22) 

2d — 1 y J 

which is located in the domain —1 < 4bA < 0. With this critical value of A, h(x,y ) 
becomes factored; 

64b 2 


(d- l) 2 ’ (d- 1) 2 

h(x, 0) = 0 also has two positive solutions x — x\(> 0) and x = x 2 (> xi), where 

—8 ak — 4b \/l + 4bA — 8ak + 4b \/l + 4bA 

X\ := -=-. x 2 :=-~-. 

3 — 4bA 3 — 4aA 

The shape of the hyperbola drastically changes at the following critical value; 

, _ r 2d - 5 

4aA. = 


h(x,y) =-(d — l) 2 ( y - 


64 kax 


+ 


x [y + 


(d — l){8b — (d — l)iu} (d — l) 2 w 
8 kwx 

+ w 


(4.23) 


8 a — (d — 1)m; 

where a constant w satisfies (d — 1 ) 2 w 2 + 16 daw + 64b 2 = 0, and therefore h(x,y ) = 0 
consists of two straight lines. 


x = x + (y) and x = x-(y) coincide when Z(y) vanishes. These points are located on 
x = xp(y) only for d = 5 or 4bA = —(2d — 5)/(2d — 1). For 4bA = —(2d — 5)/(2d — 1), 
Z(y) vanishes at 

y = -jzj=-y«- ( 4 ’ 24 ) 

For d — 5, it vanishes at 


V = 


b 20bA + 1 ± \ (4bA - 3)(36bA + 5) 


2(1 + 4bA) 

The (x, y)-planes for —1 < 4bA < 0 are drawn in Figs. W|6 


=: y 5 (±)- 


(4-25) 
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4.3.2 Existence of solutions 


In the present case, existence of static solutions depends on the value of the mass parameter 
m. First of all, as in Einstein gravity, it is sh own that there is no static solution for 
sufficiently small \m\. For —1 < 4bA < 0, Eq. (4.15) shows that x + (y ) < x mSLX (y) holds 


in the domain y < y c (+)- As seen in Fig. [2j the curve y = m/x <4 -5 )/ 2 moves to the right 
as m(< 0) increases approaching the x-axis in the limit of m —> —0. Therefore, there 
exists a critical value m c such that, for m c < m(< 0), the intersection of x — x+(y) with 
y = m/x is located outside the physical domain in the (x,y) plane and hence there 
is no static solution. This critical value m c is obtained by solving the following algebraic 
equations: 


2 x = (d — 1 )ky(m c ) — Aka , 

(d — l) 2 K.y(m c ) = 4 (d — l)bA + d — 3 + \J 4 (d — 1 )(d — 5)bA + (d — 3) 2 , 


(4.26) 


where y(m c ) := m c /x' d 5 ^ 2 . 


In the case of — (2d — 5)/(2d — 1) < 4dA < 0, Z(y) is non-negative and hence x = x + (y) 
is continuous. Therefore, as seen in Figs. [|and[5j there is a static solution for each value 
of m satisfying m < m c . 


In the case of —1 < 4aA < —(2d — 5)/(2d — 1), in contrast, Z(y) is negative in the 
domain of < y < y + , where 


y± ■= 


8a AdaA + (d — 4) ± y (4<5A — 3) [4(2d — l)aA + 2d — 5] 
(d — 1) 2 (1 + 4dA) 


(< 0). (4.27) 


As seen in Fig. [6j the curve x = x + (y) is no more continuous and does not exist in the 
domain of < y < y + . Since the curve y = m/x^ d ~ 5 ^ 2 moves to the left as m decreases, 
there exists a range of negative m G (m_, m + ) such that the curve y = m/x^ d ~ 5 ^ 2 does not 
intersect with the hyperbola h(x, y) = 0 and hence there is no static solution. This shows 
a sharp difference from the general relativistic case. 


4.3.3 Stability of solutions for d — 5 


Now let us study the stability of solutions. For this purpose, we use the following quantity: 


x + (y) -x P (y) 


S(y) — 2k {(d — 1 )(d — 3 )y - 16a} a/ Z(y) 
Ak(3 — 4dA) {(d — l)(d — 3 )y — 16b} 


(4.28) 
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Figure 3: The (x,y) plane with A = 0, a = 1, and d — 6. Thick solid curves consists of 
x = x + (y ) and x = X-(y), while thin solid curves are x = Xp(y) corresponding to P = 0. 
The dashed lines consist of x — x max (y) and x = x m i n (y). The shadowed region corresponds 
to P > 0 in the physical region in the domain of x > 0, y < 0. 



Figure 4: The (x,y) plane for —(2 d — 5)/(2d — 1) < 4<5A < 0. If 4dA is close to 
— (2d — 5)/(2d — 1), a part of x = X-(y) enters the physical region but the corresponding 
solutions are unstable. 
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X 



Figure 5: The (x,y) plane for 4aA = — (2d — 5)/(2d — 1). The intersection between 
x = x + (y) and x = x_(y) is located on x = xp(y) (thin solid curve) for any d. 


x 



Figure 6: The (x,y) plane for —1 < 4dA < —(2d — 5)/(2d — 1). At y — y±, Z = 0 (and 
hence x + = xJ) are satisfied. It is noted that the points (x + ,y + ) and (x + ,yJ) are located 
on x = Xp(y) (thin solid curve) only for d — 5. 
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where 


S(y) \={d — 3 )(d — 1) 2 (1 + 4aA)y 2 — 8(d — l)a(4daA + d — 4 )y + 128a 2 (l + 4aA) 
=4aA{(d — 3)(d — 1 ) 2 y 2 — 8 a(d — 1 )dy + 128a 2 } 

+ (d-3)(d- 1) V - 8a(d - 1 )(d- 4 )y + 128a 2 . (4.29) 

If S(y) is positive in some negative domain of y, then x + (y) > Xp(y) holds there, which 
means that P > 0 is satisfied at x = x + (y). Therefore, if there are intersections of 
x = x + (y) with y = m/x ( d-5 )/ 2 in the physical domain with S(y) > 0, the corresponding 
static solutions are stable. Since the stability of static solutions is different between d = 5 
and d > 6, we treat the case of d = 5 here and the case of d > 6 will be treated separately. 


In the case of d — 5, the solution cannot be unstable, shown as follows. For d — 5, 


S(y) = 2Z(y) is satisfied and so Eq. (4.28) becomes quite simple. Since Z(y) > 0 is 
satisfied on the curve x = x + (y), we have x + (y) > xp(y) there and the corresponding static 
solutions are stable or marginally stable. 


Marginally stable static solutions are realized at y = y 5 (±) satisfying Z(y 5 ^) = 0, where 
y$(±) is given by Eq. (4.25). Because the reality of y 5 (±) requires 4dA < —5/9, static 
solutions for —5/9 < 4aA < 0 are all stable. (See Fig. kl) 


On the other hand, for — 1 < 4aA < —5/9, the static solutions with m = y§(±) are 
marginally stable, while the solutions with m < ys(-) or y 5 ( + ) < m are stable. For y 5 (_) < 
m < y 5 ( + ), there is no solution. Figure]^ shows the case of 4aA = —5/9, in which y 0 
becomes y 5 ( + ) = y 5 (_) for d = 5. Figure bis hows the case of —1 < 4bA < —5/9, in which 
y + and y_ become y 5 ( + ) and y 5 (_) for d — 5, respectively. 


4.3.4 Stability of solutions for d > 6 


In order to discuss the stability for d > 6, we evaluate the function h(x, y) on the marginally 
stable curve x = xp (y): 


where 


h(x P (y),y) 


_ w(y) _ 

16{(d — l)(d — 3)y — 16d} 2 ’ 


(4.30) 


W (y) :=4aA {(d — 3)(d — l) 2 y 2 — 8 d(d — 1 )ay + 128a 2 } 2 

+ (d - 3) 2 {d - l) 4 y 4 - 16(d - 3)(d 2 -5 d+ 12){d - 1 ) 2 dy 3 
+ 64(d - l)(d 3 + 3d 2 - 52 d + 112)d 2 y 2 
- 2048(d 2 -d- 8)d 3 y + 16384A 4 . 


(4.31) 






We will show that the sign of h{xp{y),y) is definite in the domain of negative y for 
— (2d — 5)/(2 d — 1) < 4bA < 0, which means that x = x + (y ) does not intersect with 
x = xp(y). Then, by continuity of the curves x = x + {y) and x = xp{y) for y < 0, the sign 
of x + (y) — xp{y ) is the same as x+(0) — xp(0) and it is actually positive; 


x+(0) - x P (0) 


~2koi{\ T 4ctA) T 4ck^/(1 T 4ctA) 
3 — 4bA 


> 0 . 


(4.32) 


Therefore, x + {y) > Xp(y) is satished for y < 0 and hence the corresponding static solutions 
are stable. (See Fig. 0) 


In order to prove the definiteness of the sign of h(xp(y),y) for —{2d — 5)/(2 d — 1) < 
4bA < 0, we use the fact that W{y) is an increasing function of A. From the following two 
expressions; 


W'WIa-o 


I- r(t/)| 4S A._M = 5 

2d— 1 


= {d — 3 ) 2 (d — 1 ) 4 y 4 — 16 {d — 3 ){d 2 — 5 d + 12) {d — 1 ) 2 by 3 
+ 64 {d - 1 )(d 3 + 3d 2 - 52 d + 112)5V 

- 2048(d 2 -d- 8)b 3 y + 16384b 4 > 0, (4.33) 

4 {(d — 1 )y + 8b} 2 {{d — 1 ) 2 {d — 3 ) 2 y 2 — 32 d{d — 3 )ay + 256b 2 } ^ 

“ 2d -1 “ ’ 

(4.34) 


it is concluded that the sign of W{y) and hence the sign of h{xp{y),y) is definite in the 
negative domain of y. 


In the case of 4bA = — {2d — 5)/{2d — 1), the solution can be marginally stable. Because 
the equality in Eq. (4.34) holds only at y = —8b/(d— 1), h{xp{y), y) = 0 is satished only at 
y = — 8a/{d — 1) and the sign of h{xp{y),y) is definite elsewhere. Therefore, x + {y) > xp{y) 
and x + (y) = xp(y) are satished at y ^ —8b/(d — 1) and y = — 8a/{d — 1), respectively. 
Namely, the static solution with a critical value of m corresponding to y = —8b/(d — 1) is 
marginally stable and solutions with other values of negative m are stable. (See Fig. [5]) 


The situation is complicated for —1 < 4bA < —(2d — 5)/(2d — 1). In this parameter 
region, the solution may be dynamically unstable which shows a sharp difference from the 
general relativistic case. 


Figure [6] shows the (x, y)-plane in this case, x = x + {y ) exists only in the domains of 
y < y~ and y + < y{< 0) because Z{y) is negative in the domain of < y < y + , where y± 


are defined by Eq. (4.27) and satisfy Z{y±) = 0. From the following expression; 


128b 2 |4(2d - l)bA + (2d - 5)) 

Z(y 0 ) = - k ^- i (< 0), 


d — 1 


(4.35) 


where yo := — 8a/{d — 1), we obtain an inequality y_ < —8 a/{d — 1) < y + . 
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For our purpose, we use the fact that S(y) is an increasing function of A. Non-negativity 
of Z gives the following inequality: 


4aA > 3 — 


4 {(d — 1 )y — 8 a} 2 
{(d — 1 ) 2 y 2 — 16 day + 64a 2 } 


This lower bound gives a lower bound of S: 

s > 32 (d - 5 )ay {64a 2 - (d - 1) V} 

^ — (d — l) 2 y 2 — 16 day + 64a 2 


(4.36) 


(4.37) 


If the right-hand side is positive in some domain of y, x + (y) > Xp(y) is satisfied there and 
hence the corresponding static solutions are stable. 


Because of the inequality y_ < —8 a/(cl — 1) < y + , static solutions corresponding to 
y < y- are dynamically stable. In contrast, the static solution with y — y + is dynamically 
unstable since we have 


%+(y+) -x P (y+) 


8 (d — 5 )ay + {64a 2 — (d — 1 ) 2 y\} 

k {3 — 4dA) {(d — 1 )(d— 3 )y + — 16d} {(d — 1 ) 2 y\ — 16 day + + 64d 2 } 

(4.38) 


Figure [6] shows that the dynamically unstable solutions are realized only very close to 
y = y + and the solutions with y + -C y(< 0) become stable. All the results obtained in the 
present paper are summarized in Table [2} 


Table 2: The existence and stability of Z 2 symmetric static thin-shell wormholes made of 
pure negative tension in the GR branch with a > 0 and 1 + 4AA > 0. ”S”, ”M”, ”U” stand 
for ’’Stable”, ’’Marginally stable”, and ’’Unstable”, respectively 




Static solutions exist? 

Stability 

k = 1 

m > 0 

Yes 

U 

m < 0 

No 

- 

k = 0 

m = 0 

A > 0: No 

- 

A < 0: Yes 

M 


No 

- 

k = -1 

m > 0 

No 

- 

m < 0 

A > 0 : No 

- 

-(2d - 5)/(2d - 1) < 4dA < 0: Yes 

S 

4bA = —(2d — 5)/(2d — 1): Yes 

S or M 

— 1 < 4bA < —(2d — 5)/(2d — 1) with d = 5: Yes 

S or M 

— 1 < 4bA < —(2d — 5)/(2d — 1) with d > 6: Yes 

S, M, or U 
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5 Summary and discussions 


In the present paper, d(> 5)-dimensional static thin-shcll wormholes with the Z 2 symmetry 
have been investigated in the spherically (k = 1), planar (k — 0), or hyperbolically (k = 
— 1) symmetric spacetime in Einstein-Gauss-Bonnet gravity. For our primary motivation 
to reveal the effect of the Gauss-Bonnet term on the static configuration and dynamical 
stability of a wormhole, we have studied the stability against linear perturbations preserving 
symmetries in the simplest set up where the thin shell is made of pure negative tension, 
which satisfies the null energy condition. 

In this system, the dynamics of the shell can be treated as a one-dimensional potential 
problem characterized by a mass parameter m in the vacuum bulk spacetime for a given 
value of d, k, the cosmological constant A, and the Gauss-Bonnet coupling constant a. We 
have studied solutions which admit the general relativistic limit a —>■ 0 and considered a 
very conservative region in the parameter space. The shape of the effective potential for the 
shell dynamics clarifies possible static configurations of a wormhole and their dynamical 
stability. 

As seen in Tables [l] and [2j the results with and without the Gauss-Bonnet term are 
similar in many cases. For k — 1, static wormholes require m > 0 and they are dynamically 
unstable. For k — 0, static wormholes require m — 0 and A < 0 and they are marginally 
stable. For k — — 1, m < 0 and A < 0 are required for static wormholes. 

We have clarified the effect of the Gauss-Bonnet term on the stability in a perturbative 
method by expanding the equation in a power series of a. We have shown that, for a/a 2 E <C 
1, the Gauss-Bonnet term tends to destabilize spherically symmetric thin-shcll wormholes 
(k = 1), while it stabilizes hyperbolically symmetric wormholes (k = —1). For planar 
symmetric wormholes (k — 0), the Gauss-Bonnet term does not affect their stability and 
they are marginally stable, same as in Einstein gravity. However, we have observed that 
the non-perturbative effect is quite non-trivial. 

Notable difference between Einstein gravity and Einstein-Gauss-Bonnet gravity appears 
in the case of k — — 1. In Einstein gravity, static wormholes exist when the mass parameter 
m is less than a critical negative value and they are dynamically stable. This is also the 
case in Einstein-Gauss-Bonnet gravity for — (2d — 5)/(2d — 1) < 4aA < 0. However, for 
4bA = —(2d — 5)/ (2d — 1), a static wormhole becomes marginally stable if m is fine-tuned. 
For —1 < 45A. < —(2d — 5)/(2d— 1), in contrast, static wormholes cease to exist for a finite 
range of m and furthermore, dynamical property of the wormhole is different for d = 5 and 
d > 6. For d — 5, static wormholes are generically stable but become marginally stable if 
m is fine-tuned. For d > 6, in addition to them, wormholes are dynamically unstable in a 
finite range of m. I 11 summary for k = —1, the Gauss-Bonnet term shrinks the parameter 
region admitting static wormholes and tends to destabilize them non-perturbatively. 
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As the effect of the Gauss-Bonnet term on the existence and stability of static wormholes 
has been revealed in the present paper, the effect of its dilaton coupling is now of great 
interest. Unfortunately in the presence of a dilaton, exact bulk solutions are not available 
to construct thin-shell wormholes. Nevertheless, this is a promising direction of research 
leading to understand the result in [9]. We hope that the result will be reported elsewhere. 
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A Derivation of the equation of motion for a thin shell 


In this appendix, we present the details how to derive the equation of motion for the shell 
(2.15) and (2.16) from the junction conditions (2.10). 


For the following vacuum bulk metric ( 2.7[ ); 

d-Srf =g flu dx tJj dx u = — f(r)dt 2 + /(r) -1 dr 2 + r 2 ~/ABdz A dz B , 


/(r) -.=k + gl lT^/l + ^ + 4bA], 


the non-vanishing components of the Levi-Civita connection are given by 

Id/ t = J_d J = _J_d/ 

tt 2 ar ’ tr 2/dr’ rr 2/ dr ’ 

1 


r r AB — — r flAB, T Br — -5 B, 

r 


r a bc = r a bc (z), 


(A.l) 

(A.2) 


(A.3) 


where f A BC is the Levi-Civita connection on the maximally symmetric base manifold. 


In this spacetime, the position of the thin shell is described by r = a{r) and t = T(r), 
where r is the proper time on the shell. The future directed unit tangent vector to the 

shell is 


^ d -d . d 

dxv dt a dr ’ 

of which normalization condition = — 1 is written as 

(A.4) 

1 = /(o)T ' 2 - m • 

(A.5) 
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where a dot denotes the differentiation with respect to r. The unit normal one-form n /t to 
the shell is given by 

rp/TU 1 = —adt + tdr , (A.6) 

which satisfies rp/iA = 0 and n /J n /i = 1. The vector n^{d/dx^) is pointing increasing 
direction of r. 


The (d — l)-dimensional induced metric h t j on the shell is given by 

ds^_ x = /iij(f)d£‘d£ J = -dr 2 + a(r) 2 p AB dz A dz B . (A.7) 

where £° = r. Non-zero components of the Levi-Civita connection ( d- 1 )r*. fc in this spacetime 
are 


(d-l)pr _ n ^ (d-l)pA — -A A (d-l)pA _p 

1 AB — aa l ABi 1 tB — °Bi 1 BC — 1 

a 


A 

BC- 


(A.8) 


From our definition of the Riemann tensor; 


/?/' — p) - 8 r M 4 - p 1 t* 

^ L 7 y/on- j/rr ^( 7 1 i/n • 1 


pM p* 


(A.9) 


the non-zero components of the Riemann tensor 77 l y fc , Ricci tensor 77, j, and Ricci scalar 
77 are computed to give 


K t btd = ad lBD, K A bcd = (k + a 2 )(S A c j BD - 5 A D ^ BC ), 

TZ T t — — {d — 2) TZ-ab — + (d — 3)(/c + d 2 ) 

77 =2(d - 2)- + (d - 2)(d - 3) (~ 2 + ^ Y 

a \ or a z / 


From these expressions, we obtain the non-zero components of 


P BrD =^(<* — 3)(d — 4)(i + 0 2 b BD, 

P^BCD =( d ~ 4)|aS + - 5)(fc + a 2 )j(<5' 4 c 7BD - C D j B c). 


(A.10) 
(A.11) 

(A.12) 


(A.13) 
(A.14) 


The extrinsic curvature of the shell is computed from the following definition: 

Kn :=(V^)efe^ 

= - - V >.ef e" , 

where e ^ := dx^/dt; 1 . Using 

e-df = tdr , e-df = ddr , efd£ l = 5 B dz B , 


(A.15) 


(A.16) 
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and Eq. (|A.5|) together with its derivative with respect to r, we obtain the non-zero com- 

(A.17) 


ponents of K l j\ 


K T r = ( a + 4 

fT 


K A u = f —5 A 


B 


where a prime denotes the derivative with respect to a. From the above expressions, we 
compute 


KijK* 


1 

fT 
1 


K=— [ d+ J - ) + 


f \ , (d - 2) ft 


j2J'2 


r 


a + — + (d — 2) ( — , 


fT 


J T r = ~ 


(d 2)(d — 3) fT , f 


a+ 2 




3a 


B 


J = 


(d — 2)(d — 2>)fT f 3 
3a \ i 




(A-18) 
(A-19) 
(A.20) 
(A.21) 
(A.22) 


Now we are ready to write down the equation of motion for the shell. Under the as¬ 
sumptions of the Z 2 symmetry and the form of S l j as 


S l j = diag (~p,p,p, ■■■ ,p) + diag(—cr, -a, -a,-- - , -a ,), 


(A.23) 


the junction conditions (2.10) give (2.15) and (2.16), where we used Eq. (A.5) in the 
following form: 


ff\ 2 3 4 = f_ # 

a ) a 2 a 2 


(A.24) 
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